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Abstract 

In the theoretical description of recent experiments with dilute Bose gases 
confined in external potentials the Gross-Pitaevskii equation plays an important 
role. Its status as an approximation for the quantum mechanical many-body 
ground state problem has recently been rigorously clarified. A summary of this 
work is presented here. 



1 Introduction 

The Gross-Pitaevskii (GP) equation is a nonlinear Schrodinger equation that was 
introduced in the early sixties as a phenomenological equation for the order 

parameter in superfiuid He4. It has come into prominence again because of recent 
experiments on Bose-Einstein condensation of dilute gases in magnetic traps. The 
paper ^ brings an up to date review of these developments. 

One of the inputs needed for the justification of the GP equation starting from the 
many body Hamiltonian is the ground state energy of a a dilute, thermodynamically 
infinite, homogeneous Bose gas. The formula for this quantity is older than the GP 
equation but it has only very recently been derived rigorously for suitable interparticle 
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potentials. See |Q and 0. The paper ^ goes one step further and derives the GP 
as a limit of the full quantum mechanical description. This derivation is summarized 
in the present contribution. 

The starting point of the investigation is the Hamiltonian for bosons that 
interact with each other via a spherically symmetric pair-potential v{\xi — Xj\) and 
are confined by an external potential V{x): 

N 

H = Y,{-^^ + ViS^)}+ Yl ^(l^^-^jD- (1-1) 
i=l l<i<j<N 

The Hamiltonian acts on symmetric wave functions in L2(R3^,rf3^x). The pair m- 
teraction v is assumed to be nonnegative and of short range, more precisely, v{r) < 
(const.) r~(^+^) as r — »• oo, for some e > 0. The potential V that represents the trap is 
continuous and ^(f) — > oo as — * oo. By shifting the energy scale we can assume 
that min^ V{x) = 0. 

Units are chosen so that h = 2m = 1, where m is the particle mass. A natural 
energy unit is given by the ground state energy hu of the one particle Hamiltonian 
— (^^/2m)A + V. The corresponding length unit, \/h/ (mcj), measures the effective 
extension of the trap. 

The ground state wave function \l/o(^i5 • • • 5 ^n) satisfies H'^q = E^^^^q, with the 
ground state energy = inf spec H. If f = 0, then \l/o(^i, • • • , a; at) = YliLi ^o{^i)y 
with $0 the ground state wave function of —A + V{x). On the other hand, if f 7^ the 
ground state \E'o(a;i, . . . , xat) is, in general, far from being a product state if is large. 
In spite of this fact, recent experiments on BE condensation are usually interpreted 
in terms of a function $*-^^(af) of a single x G M^, the solution of the Gross- Pitaevskii 
equation 

(-A + 1/ + 87ra|<l>n$ = A$, (1.2) 
together with the normalization condition 

[ \<^{x)\^ = N. (1.3) 

Here a is the scattering length of the potential v: 

a = hm (r - ^] (1.4) 

where Uq satisfies the zero energy scattering equation. 



—u"{r) + |t>(r)'u(r) = 0, 



(1.5) 



and Mo(0) = 0. (The factor 1/2 is due to the reduced mass of the two-body problem.) 
The GP equation ( |1.2|) is the variational equation for the minimization of the GP 

energy functional 

£GP[^] = J (|V$|' + F|$|' + 47ra|<l>|^)d3x (1.6) 

with the subsidiary condition ( |1.3|) . The corresponding energy is 

E^^{N,a)= inf £^^[^] = £^^[^^% (1.7) 

J |<I>|2=Af 

with a unique, positive The eigenvalue, A, in ( |1.2| ) is related to E^^ by 

A = dE^^{N, a)/dN = E^^{N, a)/N + Anap, (1.8) 

where 

p=^l ix)\'d'x (1.9) 

is the mean density. The minimizer ^'^^ of (|1.6|) with the condition ( |1.3|) depends on 
N and a, of course, and when this is important we denote it by 

Mathematically, the GP equation is quite similar to the Thomas-Fermi-von Weizs- 
acker equation |^ and its basic properties can be established by similar means. See 
0, Sect. 2 and Appendix A. 

The idea is now that for dilute gases one should have 

^GP^^QM pQ^(f)^|$GP(f)|2 = pGP(^)^ (1_10) 

where the quantum mechanical particle density in the ground state is defined by 



p^^{x)=N j \'^q{x,X2,. . . ,XN)\'^dX2 - ■ -dxN- (1.11) 
Dilute means here that 

pa^ < 1. (1.12) 

The task is to make ( |l.l(j| ) precise and prove it! 
The first remark is that by scaling 

E«P(Ar,a) = ArE^P(l,Ara) and = A^i/^^gp^. (1.13) 

Hence Na is the natural parameter in GP theory. It should also be noted that E^^ 
depends on A^ and v and not only on the scattering length a. However, in the limit 
we are about to define, it is really only a that matters. To bring this out we write 

v{r) = (ai/a)^t>i(air/a), (1-14) 

where Vi has scattering length ai and regard vi as fixed. Then E'^^ is a function 
E^^{N, a) of A^ and a and we can state our main result. 



Theorem 1.1 (Dilute limit of the QM ground state energy and density). 



lim ' = 1 1.15 

and 

lim ^P^!^(:?) = |$^L(^)r (1-16) 
if Na is fixed. The convergence in is in the weak Li-sense. 

In particular, the limits depend only on the scattering length oi = Na of Vi and not 
on details of the potential. 

Remark. Since p ^ N the attribute 'dilute' may seem a bit strange for these limits. 
However, what matters is that the scattering length a is small compared to the mean 
particle distance, ~ p^^^, and if Na is kept constant, then pa^ ~ A^~^. It is also 
important to remember that the unit of length, ^Jh/ {muo) = fltrap, depends on the 
external potential, and a really stands for a/otrap- If the external potential is scaled 
with N, the both Otrap and the energy unit fhuj depend on N. For instance, in order to 
achieve a finite transition temperature for the BE condensation of a noninteracting 
gas in a parabolic trap it is necessary to keep Nu^ fixed in the thermodynamic limit, 
and hence atrap ~ A^^''^ and fku ~ A^~^/^. (See 0], Eq. (14).) However, since the 
energy unit cancels in ( |1.15| ) a dependence of on does not affect the validity of 



( |1.15|) , and (|1.16|) also remains valid taking into account that both sides really contain 



the factor cttrap- 



2 The dilute homogeneous Bose gas 

The motivation for the last term in the GP energy functional ( |1.6| ) is an asymptotic 
formula for the quantum mechanical ground state energy Eo{N,L) of bosons in 
a rectangular box of side length L (i.e., the homogeneous case), that was put for- 
ward by several authors many decades ago. Consider the energy per particle in the 
thermodynamic limit with p = N/ fixed: 

eo{p)= lim Eo{pL\L)/{pL^). (2.1) 

L— >oo 

According to the pioneering work of Bogoliubov the leading term in the low density 
asymptotics of eo(p) is given by 

eo(p) ^ iupa (2.2) 



for pa^ ^1. In the 50's and early 60's several derivations of this formula were 
presented [jT^.They all depended on some special assumptions about the ground state 



that have never been proved or the selection of some special terms from a perturbation 
expansion that most likely diverges. The only rigorous estimates of this period were 
obtained by Dyson [|lT] for hard spheres: 



1 + 2^1/3 eo(p) 1 , , 

> Jlin. > (2.3) 



(1 - Fi/3)2 - 4T,pa - lOv^ 

with Y = Airpa^/S. The upper bound has the right asymptotic form, and it is not 
very difficult to generalize it to other potentials than hard spheres. The lower bound 
on the other hand is too low by a factor of about 1/14 and a remedy for this situation 
was obtained only 40 years after Dyson's paper: 

Theorem 2.1 (Lower bound for a homogeneous gas). If v is nonnegative and 
of finite range, then 

^>(1-CY'/'') (2.4) 
An pa 

with some constant C . 



The proof of this theorem is given in see also the exposition in [^. For the 
application to the proof of Theorem 1.1 we need a version for finite boxes that is 
implicitly contained in |^ and explicitly stated in 0: 

Theorem 2.2 (Lower bound in a finite box). For a positive potential v with fi- 
nite range there is a 6 > such that 

Eo{N,L)/N >A7rppa{l-CY^^^^) (2.5) 

for all N and L with Y < 6 and L/a > C'Y~^I^\ Here C and C are constants, 
independent of N and L. (Note that the condition on L/a requires in particular that 
N must he large enough, N > (const. )y~^/^''.y) 

These results are stated for interactions of finite range. An extension to potentials 

V of infinite range decreasing faster than 1/r^ at infinity is obtained by approximating 

V by finite range potentials, controlling the change of the scattering length as the cut- 
off is removed. In this case the estimate holds also, but possibly with an exponent 
different from 1/17 and a different constant. It should be noted, however, that the 
form of the error term in ( f^.4| ) is dictated by the method of proof. The true error 
term presumably does not have a negative sign for sufficiently small Y. 



3 An upper bound for the QM energy 



We now turn to the inhomogeneous gas. In order to prove Eq. ( 1.15|) one has to 
establish upper and lower bounds for in terms of E'^^ with errors that vanish 



in the limit considered. As usual, the upper bound is easier. It is based on test wave 
functions of the form 

N 

...,xn) = 1[ ...,xm). (3.1) 

i=l 

where F is constructed in the following way: 

N 

F(xi, . . . , xat) = JJ /(ti(xi, . . . , Xi)), (3.2) 
1=1 

where ti = min{|xj — Xj\,l < j < i — 1} is the distance of Xi to its nearest neighbor 
among the points Xi, . . . , and / is a function of t > 0. With uq the zero energy 
scattering solution and /o(r) = Uo{r)/r the function / can be taken as 

f{r) = (3.3) 

for r < b = (47rp/3)~^/^ and 1 otherwise. The function (|3.1| ) is not totally symmetric, 
but for an upper bound it is nevertheless an acceptable test wave function since the 
bosonic ground state energy is equal to the absolute ground state energy. 
The result of a somewhat lengthy computation is the upper bound 

<E^P(l + 0(F^/=^)) (3.4) 

with Y = Ana^p/S. Note that Y ~ N'"^ since a ~ N'^ and N. 

4 The lower bound 

To obtain a lower bound for the QM energy in an external potential the strategy 
is to divide space into boxes and use the estimate (|2.5|) for a homogeneous gas in 
each box with Neumann boundary conditions. One then minimizes over all possible 
divisions of the particles among the different boxes. This gives a lower bound to the 
energy because discontinuous wave functions for the quadratic form defined by the 
Hamiltonian are now allowed. Finally, one lets the box size tend to zero. However, it 
is not possible to simply approximate V by a constant potential in each box. To see 
this consider the case of noninteracting particles, i.e., f = and hence a = 0. Here 
^QM _ ]\jfy^^ \yy^j^ 'ualve' box method gives only as lower bound, since it clearly 
pays to put all the particles with a constant wave function in the box with the lowest 
value of V . 

For this reason we start by separating out the GP wave function in each variable 
and write a general wave function \E' as 

N 

^(fi, . . . , f;v) = n ^''''{Si)F{xi, ...,xn) (4.1) 

i=l 



This defines F for a given \1/ because is everywliere strictly positive, being tlie 
ground state of the operator —A + + 87ra|$'^^p. We now compute the expectation 
value of H in the state using partial integration and the variational equation ( |1.2|) 
for The result is 



(VI/, VI/) 

with 



^ /nLp^^'l^.) I |V.Fp+^t;(|x,-x,|)|Fp-8vrapGP(f^)|Fp 



4 = 1 



/nf=iP^^(^;^)I^P 

(4.3) 



Compared to the expression for the energy involving v]/ itself we have thus obtained 
the following replacements: 

TV N 

F(f) -87rap^^(x), and JJ cixi JJ p'^^(fi)cifi. (4.4) 

i=l i=l 

(Recall that p^^{x) = |$^^(x)p.) We have to show that the normalized quadratic 
form Q is bounded below by —inapN, up to small errors, and we use the box method 
on this problem. 

Labeling the boxes by an index a we have 

infQ(F) > inf EinfQ„(F«), (4.5) 

a 

where Qa is defined by the same formula as Q but with the integrations limited to 
the box a, is a wave function with particle number Ua, and the infimum is taken 
over all distributions of the particles with ^ Uq, = A^. We now want to use (|2.5|) and 
therefore we must approximate p^^ by a constant in each box. Let Pa,raax and Pa^minj 
respectively, denote the maximal and minimal values of p^^ in box a. We obtain 

(4.6) 

Pa, max 

where L is the side length of each box and EQ{n,L) is the ground state energy of 
n bosons in a box without external potential. Now by Theorem 1.1 EQ{na,L) ^ 
Anan'^/L^, and 

inf V(47m^/L3 - 87rp„,^axn„) > -47r Vp^ .^L^ ^ -An I (p^^)' = -47rpiV, 

(4.7) 



which looks promising, were it not for the factor Pa,min/ Pa.max in ( |4.6| ) which we would 
like to replace by 1 for small boxes. The problem is that for any fixed size of boxes 
Pa,mm/ Pa,max tcuds rapidly to zcro for boxes far from the origin. 

This problem can be solved by enclosing the whole system in a big box of side 
length R. If is the solution of the GP equation with Neumann conditions on 
the boundary of (this is the minimizer of a GP functional where the integration 
is restricted to A^), then = is bounded from below away from zero in 

A/j. Replacing in ( [4. 3D everywhere p"^^ by and restricting the integrations to 
An we can let the side lengths of the small boxes tend to zero and be sure that 
Pa,Tam/ Pa,ma.x 1 Uniformly for all the small boxes. 

For this method to work, however, we must control the error made by enclosing 
the system in the big box. Let E^^{N,a) denote the quantum mechanical ground 
state energy in A/j with Neumann conditions on the boundary. The essential step is 

Lemma 4.1. There is an Rq < oo, depending only on Na such that 

E'^^{N,a)> Ef^{N,a) (4.8) 
for all R> Rq and all N, a with Na fixed. 

This lemma follows from V{x) — > oo for oo, together with an estimate for 

the chemical potential from above: 

Ef'iN + 1, a) - Ef{N, a) < e{Na){l + 0(F^/^)) (4.9) 

where e{Na) depends only on Na and is independent of R. The proof of ( [4.9|) is 
similar to the proof of the upper bound ( |3.4| ) but this time one uses 

<(a;i,...,f^)/(t^+i) (4.10) 

with the ground state wave function for A^-particles and / as in (|3.3|) as a test 
wave function for the + 1-particle Hamiltonian. Besides the lemma one needs 
the easy to verify fact that the GP energy E'^^ in the box A^ converges to E*^^ as 
i? — i> oo. 

The lower bound, and hence the main result, now follows from a rigorous version 

1/17 

of (|4.7| ) in the box A/j. We use (|2.5|) in the boxes a; the error here is 0{Ya ) with 



Ya ~ a^ria/L^ < a^N/L\ For Na fixed, this is 0(A^-2/i72^-3/i7)_ rj.^^ difference 
between Pa,ma.x and Pa,mm introduces an additional error, 0{L). Both error terms are 
of the same order for — > oo if we choose L ~ N~^^^^. We thus obtain 

E'^^{N,a) > Ef{N,a) > E'i^{N,a){l - {const. )N~^/^^). (4.11) 

The constant may depend on R, but this is of no harm, since we first take — > oo 
and then R —>■ oo, using E*^^ —>■ E'^^ and ( |1.13| ). This proves ( |1.15| ). Convergence of 



the densities, ( |1.16|) , is obtained in a standard way (see, e.g., |Q) by variation with 



respect to the external potential. 



5 Concluding remarks 



Bose Einstein condensation in the ground state is a concept that involves the full one 
particle density matrix 



N 



j ^'o(^, ^2, • • • , Xn)*'^o{x\ X2, . . . , XN)dX2 ■ ■ ■ dx^ 



(5.1) 



and not only its diagonal, p^^{x) = 7Ar(x, x) considered here. Expressed in terms of 
creation and annihilation operators we can also write 



7Ar(x, x') = (a*(x)a(x'))o, 



(5.2) 



with {■)o the expectation value in the ground state. As an integral kernel, is a 
positive trace class operator. If Nq denotes its largest eigenvalue BE condensation 
means that there is a constant c > so that 

A^o > cN (5.3) 

for all A^. This definition applies also if is replaced by a thermal equilibrium state 

at nonzero temperature. 

Note that the density in momentum space is 



a*{p)a{p))o = / exp{ip- {x — x))'y]\f{x^ x')dxdx' 



(5.4) 



and this differs from 



/ exp{ip- x)^y p{x)dx 



unless 7Ar has rank 1. 



It is often claimed that ^'^^ is (approximately) the eigenfunction to the highest 
eigenvalue of 7Ar(x, x') and hence that |$*^^(p)p gives the momentum distribution of 
the condensate, but this is not proved yet. In fact, so far the only cases with genuine 
interaction where BE condensation has been rigorously established in the ground state 
are lattice gases at precisely half filling. (The hard core lattice Bose gas corresponds 
to the XY spin 1/2 model and BE condensation was proved in dimension > 3 |jl2| 



and dimension 2 [T^. The hard core lattice gas with nearest neighbor repulsion 
corresponds to the Heisenberg antiferromagnet and condensation was proved for high 



dimension in |12] and dimension > 3 in [14]. Dimension 2 is still open. 
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